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Abstract 

We study the regime of anticipated synchronization in unidirectionally coupled model neurons 
subject to a common external aperiodic forcing that makes their behavior unpredictable. We show 
numerically and by implementation in analog hardware electronic circuits that, under appropriate 
coupling conditions, the pulses fired by the slave neuron anticipate (i.e. predict) the pulses fired 
by the master neuron. This anticipated synchronization occurs even when the common external 
forcing is white noise. 

PACS numbers: 05.45.-a, 05.40. Ca, 42.65.Pc, 42.65. Sf 
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Synchronization of nonlinear systems is a fascinating subject that has been extensively 
studied on a large variety of physical and biological systems [|I]]. While synchronization of 
oscillators goes back to the work by Huygens, the last decade has witnessed an increased 
interest in the topic of synchronization of chaotic systems 

Recently, Voss [Q] has discovered a new scheme of synchronization, called "anticipated 
synchronization". Voss has shown that by using appropriate delay lines it is possible to 
synchronize two unidirectionally coupled systems in such a way that the slave system, y{t), 
predicts the behavior of the master system, x(t). Two coupling schemes were considered: 
complete replacement, 

x(t) = -ax(t) + f (x(t - r)) 



and delay coupling, 



(1) 

y{t) = -ay(t)+f(x(t)) 



m = f(x(t)) 

m = f(y(t)) + K[x(t)-y(t-r)]. 



f (x) is a function which defines the autonomous dynamical system under consideration, K is 
the coupling strength and r is a delay time. It is easy to see that in both schemes the manifold 
y(t) = x(t + r) is a solution of the equations, and Voss has shown that in both schemes this 
solution can be structurally stable. This is more remarkable when the dynamics of the 
master system x is "intrinsically unpredictable" , as it is the case of a chaotic system. While 
in the scheme of complete replacement the anticipation time r can be arbitrarily large, the 
delay coupling scheme requires some constrains on r and K for the synchronization solution 
to be stable |Q. Despite this fact, the delay coupling scheme is more interesting since the 
anticipation time r can be varied without altering the dynamics of the master system x. 

In this Letter we study numerically and experimentally the regime of anticipated synchro- 
nization in excitable non- autonomous systems. In our case the intrinsic unpredictability of 
the behavior of the dynamical system x does not arise from a chaotic dynamics, but rather 
from the existence of an external forcing with some element of randomness. We consider 
the coupled systems 

x(t) = f(x(t))+I(t) 

m = f(y(t)) + I(t) + K[x(t)-y(t-r)], 

where l{t) represents a common external forcing. Notice that y{t) = x(i(: + r) is no longer 
an exact solution of the equations [except in the particular case of a periodic forcing I(t + 



r) = I(t)]. We show that under appropriate couphng conditions there can be a very good 
correlation between y(t) and x(t + r) which, in practice, allows the prediction of the future 
behavior of x()f:) with a high degree of accuracy. Of course, this result is more remarkable 
when the external forcing is a random signal. 

Specifically, we have considered models of sensory neurons. Sensory neurons transform 
external stimuli signals as pressure, temperature, electric pulses, etc., into trains of action 
potentials, usually referred to as 'spikes' or 'firings'. Their behavior is typical of excitable 
systems: if the forcing is above a certain threshold, the neuron fires a pulse, and after the 
firing, the recovery process produces an absolute refractory time during which a second firing 
cannot occur. 

In general, sensory neurons work in a noisy environment. As a consequence, the time 
intervals between spikes contain a significant random component, and random spikes often 
occur even in the absence of stimuli. The topics of synchronous oscillations and noise have 
received much attention (see, e.g., since it has been suggested that synchronous firing 
activity of sensory neurons might be a part of higher brain functions and a method for 
integrating distributed information into a global picture 0. 

Here we find that the interplay of coupling, delayed feedback, and common noise can lead 
to anticipated synchronization. We illustrate this effect in the well known FitzHugh-Nagumo 
and Hodkey-Huxley neuron models. By coupling two of such systems in an unidirectional 
configuration as in the scheme (^), we find that when both systems are subjected to the 
same external random forcing, the slave system fires the same train of spikes as the master 
system, but at a certain amount of time earlier, i.e., the slave predicts the response of the 
master. 

First we show results based on the FitzHugh-Nagumo model. It consists of two variables 
x = (xi,X2). The fast variable, Xi, is associated with the activator, and the slow recovery 
variable, X2, is associated with the inhibitor. The equations for the master {xi,X2) and the 
slave y = {yi, 2/2) systems, under unidirectional coupling are, respectively (see the schematic 
diagram shown in Fig. |l]): 

Xi = -Xi{xi- a){xi-l) - X2 + 1(1) 

(4) 

X2 = e{xi - bx2) 
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and 

m = -yiivi- a){yi-l) -y2 + I{t) + 

+K[xi{t) - yi{t - r)] (5) 
?/2 = e(?/i - by2) 

where a, b, and e are constants, K is the couphng strength and r is a delay time (associated 
to an inhibitory feedback loop in the slave neuron). Note that only the fast variables of 
the two systems are coupled. When the common external forcing, I{t), is constant in time, 
the anticipated synchronization manifold: Xi{t + r) = yi{t), X2{t + r) = y2{t) is an exact 
solution of Eqs. (^) and (^). If the external forcing is above threshold and for appropriate 
values of K and r, the master system fires pulses periodically and the coupling induces a 
constant time shift r between master and slave spikes. 

We have considered different types of random external forcing I{t). The first one cor- 
responds to a random process whose amplitude remains constant for a time T and then it 
switches to a new random value chosen uniformly in [Iq — D,Iq + D], where D is the noise 
intensity. We chose Iq very close to (but below) the firing threshold of the excitable system. 
It would appear at first thought that with this type of external forcing the behavior of the 
master system can be easily predictable. However, there are two main factors that make 
the system response unpredictable: if the effect of the perturbation is not strong enough 
the system does not fire a pulse; besides, the system has a refractory time (after firing a 
pulse) during which, another firing is not possible. Figure |^ shows that anticipation occurs 
with this type of random external forcing for an appropriate value of the coupling strength 
K: after an initial transient time the two systems synchronize such that the slave system 
anticipates the fires of the master system by a time interval r. The firings in the master 
and the slave systems start at about the same time, and the anticipation phenomenon grows 
during the rising of the pulse. When the master system noisily evolves near the stable point, 
the anticipation vanishes. In other words, anticipation is a local process, during firings. 

The same qualitative results are found with other types of external forcing such as colored 
or even white noise. Figures ^(a-b) display the spikes of the master and slave systems when 
I{t) is Gaussian white noise. 

Sometimes the slave system makes an error in anticipating the master firings. While the 
slave system always fires a pulse when the master system fires a pulse, it also might fire a 
"extra" pulse, which has no corresponding pulse in the train of pulses fired by the master. 
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Notice that in Fig. |^(a) an error at about t = 1900 occurs. Not surprisingly, we find that 
the longer the anticipation time r, the larger the number of errors. However, for a given 
anticipation time, the number of errors can be reduced considerably if a "cascade" of an 
adequate number of slave neurons is considered. A detailed study of the number of errors 
and its dependency with the type of external forcing will be reported elsewhere. 

Next we show simulations based on a more realistic model, namely the model of electro- 
receptors proposed by Braun et. al This model is a modification of the Hodgkin-Huxley 
neuron model: Cmx = —iNa — — isd — isr — k, where x is the potential voltage across 
the membrane and Cm is the capacitance; iNa and ix are the fast sodium and potassium 
currents, igd and igr are additional slow currents, ii is a passive leak current. For details and 
functional dependence of the currents on the voltage x and other factors (as temperature) 
see i. 

We extend the model to account for two unidirectionally coupled neurons, with a delayed 
feedback loop in the slave neuron, and subject to a common external forcing I{t), in the 
same way as in the FitzHugh-Nagumo model, e.g., the equations for the master, x, and for 
the slave, y, neurons are: 

CmX = —i^Na ~ "^K ~ '^sd ~ ^sr "~ + H^) 

CmV = -lla - - ^Id - ^Ir - + I{t) (6) 

+ K[x{t) - y{t - r)] . 

Figures 0(c-d) display the results when the common external forcing I{t) is a Gaussian white 
noise. We chose parameters such that in the absence of forcing there are no spikes (sub- 
threshold, noise-activated firing regime). The behavior observed is qualitatively the same as 
in the FitzHugh-Nagumo model (the slave neuron anticipates the fires of the master neu- 
ron), which indicates that the anticipation phenomenon is general and model independent. 
Remarkably, in this model the anticipation time can be larger than the pulse duration. It 
is worth mentioning that anticipated synchronization is also observed in this model for pa- 
rameters such that there is spontaneous (regular or irregular) spike activity (suprathreshold 
firing regime). 

To assess the robustness of the anticipated synchronization observed in the numerical 
simulations, we have implemented the FitzHugh-Nagumo model in analog hardware and 
constructed two coupled electronic neurons (a simplified version of the circuit is shown in 
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Fig. ^). The electronic neurons were built using operational amplifiers and the cubic non- 
linearity described by x{x — a){x — l) was implemented using analog multipliers (AD633) in a 
circuit not shown for simplicity. The resistor Rc controls the strength of the unidirectional 
coupling between the master and the slave neurons. The resistor R^, [Ru = Rc in our 
case) controls the strength of the delayed feedback into the slave neuron. The coupling and 
the delayed feedback have opposite signs: while the master signal was obtained at point 
B of Fig. ^, where the voltage is —Vm, the slave signal that goes into the delay line was 
obtained at point C, where the voltage is +K. The different signs are due to the inverters 
that are located in between points A and B and C and D. The threshold on both neurons 
was controlled by a potentiometer represented by its equivalent circuit: offset and Rq. The 
analog delay line for the delayed feedback in the slave neuron was built using bucket brigade 
circuits (MN3004). A function generator with white noise output capabilities (HP33120A) 
was used to excite both electronic neurons. The signals were acquired using Lab View and 
National Instruments DAQ 6025E data acquisition board. 

Similar electronic neurons have been implemented in 0, where it was shown that their 
behavior is very similar to that of biological neurons: when interfaced to biological neurons, 
hybrid circuits, with the electronic neurons taking the place of missing or damaged biological 
neurons, could function normally. 

Our electronic coupled neurons behave very similar as in the numerical simulations. For 
an appropriate value of the coupling resistance Rc, we observe that, after a transient, 
the master and slave electronic neurons synchronize in such a way that the slave neuron 
anticipates the fires of the master neuron by a time interval approximately equal to the 
delay time r of the feedback mechanism. Figure |^ (a) shows a typical spike train, and Fig. 
^(b) displays in detail a single spike [§. We observe that, as in the numerical simulations, 
the firings of the master and the slave neurons start at about the same time: anticipation 
begins during the rising of the peak and it vanishes when the neurons are in the unexcited 
state. Without coupling and feedback {Rc = Rd = 0) the neurons fire pulses which are, in 
general, desynchronized (due to the small mismatch between the circuits). 

To summarize, we have studied the regime of anticipated synchronization in coupled sys- 
tems exhibiting neuronal-type excitable behavior, when they are driven by common external 
aperiodic forcing. We have shown that under appropriate conditions, the slave system can 
anticipate the random spikes of the master system. This is despite of the fact that the 
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anticipated synchronization manifold is not a solution of the equations. We have simulated 
numerically the coupled neurons with the FitzHugh-Nagumo and a modified Hodgkin-Huxley 
models and we have considered different types of random forcing. The FitzHugh-Nagumo 
model was also implemented in analog hardware, showing that the anticipation phenomenon 
is very general and robust. 

Our results show that non-linearity, noise and delayed feedback might conspire to produce 
new interesting and unexpected phenomena, and we hope that our findings will stimulate 
the search for anticipated synchronization in biological systems. 

The work is supported by MCyT (Spain) and FEDER, projects BFM2001-0341-C02-01 
and BMF2000-1108. C. MasoUer acknowledges partial support from the UIB. 
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FIG. 1: Schematic diagram of two model neurons coupled in a unidirectional configuration, sub- 
jected to the same external forcing and with a feedback loop (with a delay time r) in the slave 
neuron. 
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FIG. 2: Anticipated synchronization from a numerical integration of the FitzHugh-Nagumo set of 
Eqs.(|-|). The parameters are: a = 0.139, b = 2.54, e = 0.008, K = 0.15. The external forcing I{t) 
(displayed in the lower panel) is a random amplitude noise of period T = 2, mean value /q = 0.03 
and amplitude D = 0.01. Notice (upper panel) that the pulse of the slave system yi{t) (dashed 
line) anticipates the pulse of the master system xi{t) (solid line) by a time approximately equal to 
the time delay r = 4. 



9 








> 
B 

^ -4 
> 

— 6 





1 1 1 1 




1 M M 


(c) 






(d) ; 


: i 

i L : 

1 f; 





i 


1 , ' , 


1 

L k 1 









1000 3000 5000 1650 1750 
time [ms] time [ms] 



FIG. 3: Trains of spikes obtained from numerical simulations of models of unidirectionally coupled 
neurons subjected to the same external forcing, which is a Gaussian white noise with mean Iq 
and correlations {[I{t) — Io][I{t') — Iq]) = 2D5{t — t'). (a) Simulation of two FitzHugh-Nagumo 
neurons, Eqs. (|-|). The parameters are a = 0.139, h = 2.54, e = 0.008, /q = 0.03, K = 0.03, 
r = 10, D = 2.45 x 10~^. (b) Simulation of two Hodgkin-Huxley neurons, Eqs. (|6|) with K = 0.03 
ms~^, T=50 ms, and L'=0.5 mV^/ms; all other parameters as in Q (T = 6 C, V; = —75 mA in 
the notation of that paper). Left panels show typical spike trains; right panels show with detail a 
single spike. The solid (dashed) line represents the output of the master (slave) system. 
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FIG. 4: Circuit implementation of two coupled neurons. Ri = 125 kU, R2=50 kJ7, i?3=10 kJ7, 
R(. = Rjj=lOO kn, Rf=10 kn, Rn=10 kn, Ro= lO kn, Ci=100 nF, C2=l /xF. 
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FIG. 5: (a) Experimental train of spikes that shows anticipation in the spikes fired by the slave 
neuron (upper trace) with respect to the spikes fired by the master neuron (lower trace), (b) 
Detail of a spike fired by the master neuron and anticipated spike fired by the slave neuron. The 
anticipation time is 14 ms approximately. 12 



